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ON THE DIURNAL TIDE WITHIN THE THERMOSPHERE 

H. Volland and H. G. Mayr 
Goddard Space Flight Center, Greenbelt, Md. 


ABSTRACT 

The atmospheric tide within the thermosphere driven by the solar EUV heat 
input has a predominant diurnal component, contrary to the tides within the 
lower atmosphere where the semidiurnal component prevails. Based on this 
observatorial fact, a theory of the thermospheric tides is given in which the 
spherical harmonics of lowest degree are considered as the eigenfunctions of 
the problem. Perturbation theory is used which leads to a complete separation 
between a zonal and a diurnal system. A general solution is given and the zonal 
and the diurnal horizontal windsystem at 300 km height is calculated. This 
windsystem has an equatorward directed component at low latitudes during day- 
time which allows to explain the equatorial F2 anomaly as resulting from an 
enhancement of electrons at ±10° latitude due to such winds. 



ON THE DIURNAL TIDE WITHIN THE THERMOSPHERE 


1 . Introduction 

While the knowledge as well as an adequate basic theory of the tidal motion 
of the lower atmosphere are more than 100 years old [see e.g. the review articles 
of Wilkes (1951), Kertz (1957) and Siebert (1961)] , the diurnal tidal variation of 
the thermospheric density has been discovered only _0 years ago fi’om satellite 
drag measurements (see e.g. the review article by Priester, et al, 1967). The 
most important difference between the atmospheric tides within lower and upper 
atmosphere is the predominance of the semidiurnal component in the lower 
atmosphere and the predominance of the diurnal component in the upper atmos- 
phere. 

There exists general agreement that the driving force within the whole 
atmosphere is predominately solar radiation — absorbed by water vapor in the 
troposphere, by ozon in the stratosphere and mesosphere and by oxygen within 
the thermosphere (Lindzen, 1967). Thus, the diurnal component of the driving 
force is much larger than the higher harmonics. The predominance of the 
semidiurnal tide within the lower atmosphere in spite of its smaller excitement 
has originally been considered as a resonance effect. One eigenvalue of the 
atmosphere comes very close to 12 hours, the period of the semidiurnal solar 
tide. Recently, however, suggested by Siebert (1961) and quantitatively worked 
out by Kato (1966) and Lindzen (1967), it has been shown that the diurnal com- 
ponent really is the extraordinary or.e being an evanescent mode with negative 
equivalent depth so that its propagation is suppressed. It will be shown in 
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section two of this paper that the diurnal solar tidal wave becomes a propagation 
mode within the thermosphere, thus predominating the tidal motion there. 

Tidal theory of the thermosphere started with the one dimensional model of 
Harris and Priester (1962) in which heating within a vertical column and its re- 
sultant dynamic behavior have been studied. The main shortcomings of this theory 
were that the calculated diurnal density variation was not in agreement with the 
observations. Harris and Priester therefoi'e introduced an ad hoc second heat 
source in order to shift amplitude and phase of the calculated values toward the 
measured data. It has been shown from a two dimensional equatorial model 
(Volland et al., 1969) that the natural response time of the thermosphere 
with respect to maximum solar heating is just 2 hours in agreement with the 
observations and that a diurnal tidal gravity wave propagating from the lower 
atmosphere into the thermosphere is superposed to the diurnal waves generated 
within the thermosphere by the EUV heat source. 

A complete theory of the diurnal tide within the thermosphere is more com- 
plicated than the classical tidal theory because heat conduction has to be taken 
into account. Therefore, the basic tidal differential equation becomes a fourth 
order equation contrary to the second order Laplace equation in the classical 
case. Moreover, the number of eigenvalues doubles, and there does not more exist 
a clear distinction between propagation modes and evanescent modes because the 
eigenvalues are complex. On the other hand, the theory of the thermospheric 
tides in some ways becomes easier to handle than the classical theory because 
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of the predominance of the diurnal component and because the boundaries of a 
given thermospheric model behave like free internal boundaries. 

Subject of this paper is to give the basic equations and general solutions of 
the thermospheric diurnal tide. Furthermore, we shall discuss the implications 
of these solutions on the thermospheric wind system and calculate the horizontal 
wind field at 300 km altitude. An exact numerical integration of the tidal motions 
will be given in a subsequent paper. 

2. Propagation Modes and Evanescent Modes 

In this section we will give a very crude description of the propagation 
modes and the evanescent modes in tidal theory based on a plane model of the 
earth and neglecting Coriolis force and heat conductivity. We do this merely to 
achieve a basic physical picture of the nature of the wave modes and its change 
within the different atmospheric levels . 

We idealize the earth's atmosphere during equinox by a rectangular trough 
of length = 2-7rR and of width \ x . R is the earth's radius, x represents the 
meridional component, positive in south direction, y represents the longitudinal 
component, positive in east direction. The length of A. x will be defined below. 

The lower boundary of the trough is the earth's surface. The trough is open in 
the upward direction. The atmosphere within this trough shall be isothermal 
and shall behave adiabatically. In such atmosphere the diurnal solar heating can 
generate gravity waves <">f the following type: 

C nm = A nm * COS ( k nn,y) * eX PO T ± J A nm Z ) W 
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with 


fn (n + 1) - m 2 


nm \ 


R 


x -component of horizontal wave number 


2 77 m m 

k ™ = T = ~R 

y 


y -component of horizontal wave number 


/n(n + 1) 


a 


2 + k 2 horizontal wave number 

m nm 


w m = m fi angular frequency of the wave mode 
Q angular frequency of the rotating earth 

r = t + y local time 

CO J 

m 

t universal time 

x, y,z coordinates in jouth-, east- and upward direction 
(m,n) integers determining the domain of wave mode. 


A 


nm 



( 2 ) 


eigenvalue of the gravity wave of domain (n,m) [see Hines (1960), Equation (14); 
and Siebert (1961), Equation (4.12)] 


C velocity of sound 


yg 

2C 


ffl b 


2 jy - 1 
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Brunt-Vaisala frequency 
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g gravitational acceleration force 

y ratio between the specific heats at constant volume and at constant pressure. 


g = 2C ^ n(n + 1) (y - 1) 

nm mR / 


( 3 ) 


The last approximation in Equation (2) is valid because within the ranges of fre- 
quency and height, which we consider, it is « a >> co . 


The plus sign of the exponential factor in Eqi aon (1) is related to an upward 
propagating gravity wave, the minus sign is related to a downward propagating 
wave. We define as propagation modes waves with real eigenvalues A nm [S nm > 1 
in Equation (2)] and as evanescent modes waves with imaginary eigenvalues 
A m = j | A nm | [S nm < 1 in Equation (2)] . Evanescent waves are suppressed very 
rapidly during their propagation. They do not participate in an effective wave 
energy flux. A wave changes its character from a propagation mode into an 
evanescent mode at S nm = 1. The boundary condition of zero wind velocity at the 
ground (z = 0) determines the remaining amplitude factors A nm of a combination 
of an upgoing and an downgoing wave of the same domain (n, m). 


The specific choice of the x -component of the horizontal wavenumber in 
Equation (1) has been made in order to allow a direct comparison with the eigen- 
values on a sphere. These are for a nonrotating spherical earth the spherical 
functions P nm ( 6 ), where P nm are the associated Legendre's polynoms and 6 is 
the colatitude. This comparison can be made by replacing the term cos (k nm y) 
in Equation (1) by P nm . In the classical theory of tides the horizontal wave- 
number k n is equal to 
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( 4 ) 


k = — 

" ^ 

where h n is the so-called equivalent denth, because water waves in shallow 
oceans of depth h n have in fact the horizontal phase velocity 



According to our definition, k n is the horizontal wavenumber which deter- 
mines the horizontal phase velocity V n of thp characteristic waves of frequency 

&L. 

m 

Table I contains the numbers S n j of the first three symmetric functions P n 1 
from Equation (3) for the diurnal modes (m = 1). Ro,v 1 in Table I has been cal- 
culated for an isothermal lower atmosphere with temperature T 0 = 250° (C = 

320 m/sec). We notice that the eigenfunction P u has a value Sj x < 1 and thus 
belongs to the evanescent modes, while the higher eigenfunctions P 31 and P S1 are 
propagation modes. In an exact calculation taking into account Coriolis force 
the eigenfunctions P nm must be replaced by the so-called Hough functions 

CD 

e = / c< n > p, 

nm / . f?m 

which are sums of the associated Legendre's polynoms centered around the pre- 
dominant function P„ m . The determination of the eigenvalues of 6 is much 
more complicated than in the case of the nonrotating earth. Here, imaginary 
horizontal wavenumbers k n and therefore negative equivalent depths h n appear 
in Equation (2) for 6 l , and (9 3l (Lindzen, 1967). Thus, the condition of an 
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evanescent mode for P lt is even stronger valid than in the case of a nonrotating 
earth. It is this fact, the evanescent mode type of the function 

6 ~ P 

w \\ '•'u r u • 

which causes its weak generation within the lower atmosphere in spite of the 
predominant exitation force which goes into this wave domain by the solar heating. 

Table II, row 2, contains the numbers S n t calculated for a temperature cf 
T 0 = 1000°K (C = 750 ir./see) which is a typical thermospheric temperature. 

Now we notice that evei the first eigenfunction P n belongs to the propagation 
modes. Heat conduction and ion drag reduce the influence of the Coriolis force 
within the thermosphere (Volland, 1969). On the other hand, the difference be- 
tween propagation modes and evanescent modes becomes vague because the 
eigenvalues are complex giving rise to wave energy dissipation even of propaga- 
tion modes. Full wave calculations of the propagation of free internal gravity 
waves within the thermosphere show that maximum transmission of gravity 
waves exists near S nm = 1 and that with increasing S nm > 1 the transmission 
coefficient of gravity waves drops like 

IT ! ~ (S V b z (S > n 

I nm! ' w nui' ' 

where b is a positive number depending on frequency and the model adopted 
(Volland, 1968b). 

We expect therefore that the high temperature, the specific geometric and 
rotational data of the earth and the influence of heat conduction and ion drag 
cause the preference of the function P t . in the diurnal tidal wave propagation 
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within the thermosphere. This is exac. what we notice from the observed 
data of density and temperature (see e.g. Jacchia and Slowey, 1967; Priester et 
al., 1967, Taeusch et al., 1968). The predominance of the term P 1X begins at 
about 110 km altitude as we can see from the analysis of the geomagnetic S q 
current (Kato, 1956). 

3. The Basic Equations 

The complete basic equations of conservation of mass, momentum and energy 
and the equation of state written in spherical coordinates are given elsewhere 
(Priester et al., 1967). The following analysis makes use of some assumptions; 

a. Application of the perturbation theory and use of a given mean 
atmospheric model 

Perturbation theory implies that all physical parameters, depending on time 
or latitude, are small compared with the mean values averaged over time and 
sphere. Thus, higher order terms of these parameters shall be neglected. This 
is a sufficient approximation as long as the relative magnitudes of the coefficients 
in the series of spherical harmonics are smaller than 0.3: 

l a nml 

< 0.3 (n > 0) , 

a oo 

This condition holds in fact for all parameters of the thermosphere below about 
400 km altitude (see e.g. CIRA, 1965; Volland et al., 1969). Numerical cal- 
culations, which will be presented elsewhere, therefore, are limited to this 
height range. 




An important result of perturbation theory is, that the variations can be 
complete^ separated into domains of numbers m , because any coupling between 
the different domains m only occurs via nonlinear terms of the perturbation vari- 
ables. The validity of perturbation theory, on the other hand, justifies the static 
diffusion model of Jacchia (1964). Therefore, we can use mean values of tem- 
perature, density and molecular weight from the Jacchia mode]. Molecular 
weight as well as the coefficient of heat conductivity show small time and lati- 
tudinal variations below 400 km altitude in the Jacchia model. These values 
therefore are considered to be only height dependent. The collision number 
between ions and neutrals is proportional to the ion density. This value strongly 
depends on height, latitude and time. In order to be consistent in our model we 
can only take into account its height and latitude dependence. 

b. Use of a given independent EUV heat source 

The EUV heating of the neutral atmosphere via inelastic collisions with the 
ions depends itself on the temporal state of density, temperature and composition 
of the neutral air. For convenience we shall use a given independent EUV heat 
source which can e.g. be derived from the calculations of Harris and Priester 
(1962) and extrapolated into the whole sphere. In view of the large uncertainties 
in the determination of the efficiency factor of the heating this does not add 
very much to the errors already involved in the problem. 
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c. Negligence of molecular viscosity 


Molecular viscosity mainly influences the horizontal wind field via vertical 
windshear. It can be shown that below 400 km its influence on the dynamics of 
the diurnal tides within the thermosphere is small compared with the accelera- 
tion force or the ion drag force in the equation of conservation of horizontal 
momentum (Geisler, 1967; Kohl and King, 1967; Volland et al., 1969). The 
negligence of viscosity greatly simplifies the handling of the differential equations. 


d. Wind components are perturbation variables 

We shall treat all wind components including a mean longitudinal wind v 00 
as perturbation variables, thus neglecting squares and higher order terms. This 
is justified because generally the condition 


holds below 400 km altitude. a nm is the coefficient of one of the wind components 
within the series of spherical harmonics, and C is the velocity of sound in the 
particular height. Furthermore, we shall assume that the ions are completely 
bound by the earth's magnetic field, but can move freely with the neutral wind 
along the geomagnetic lines of force. Then the relation between ion velocity v i 
and neutral air velocity v is given by 

(v • B 0 ) B 0 

v. = (5) 

where B 0 is the magnetic field of a dipole which approximates the geomagnetic 
field. 
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e. Predominance of lowest domain numbers n 


Our model is valid for equinox conditions. Therefore, the asymmetric 
spherical harmonics of density and temperature variation disappear. Moreover, 
w-‘j confine ourselves to the diurnal component of the tides. We treat, therefore, 
only the domains u - 0 and m = 1 of number m. As outlined in section 2 we have 
reason to believe from a theoretical point of view that the term P x x is predomi- 
nant throughout the theimosphere. This is confirmed by the measurements 
(Priester et al., 1987; Jaechia and Slowey, 1967; Taeusch et al., 1963). Likewise, 
the density observations show a predominance of the term P 20 in the meridional 
variation (Newton, 1968). We, therefore, assume that within the domains 
m = 0 and m = 1 the higher order n~ terms of the density variations are small 
compared with the perturbation coefficients of P 20 and Pj t and can be neglected. 
Equivalent assumptions are made for the other physical parameters. This 
classification of the coefficients of the spherical harmonics into three classes 
of successive importance: 


a oo » 


! a 20 I >> l a nol 


(n A 4) 


a 00 >> Kli >:> l 3 nll ( n ^ 3 > 


( 6 ) 


greatly simplifies the calculations. It means that we can in fact treat the spheri- 
cal harmonics P 20 and P lt as the eigenfunctions of systems m = 0 and 1, re- 
spectively. 

With the assumptions outlined above the system of equations governing the 
dynamic behavior of the thermosphere is given by 
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3 u 

Tt 


+ v sin 2 I u - 2fl cos 9 v ~ v sin I cos I 


w } 


3p 
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( 8 ) 


oo "( 


3 v 
Poo 1 Tt 


+ 2fi cos P u + v v + 2Q sin 0 w r + 


} 


r sin 6 3 X. 


{ 9 ) 


/ 3w 
Poo |_3t 


v sin I cos I u - 2fi sin 0 v + v cos 2 I w V + g,o + 


} 


d p _ 
3 r 


= 0 ( 10 ) 


3 T ^ T oo 1 3 / 3 T\ 

Poo c v Tt + p oo c v w "TT" ~ 7 2 37 ( r *00 97 j ' *00 AT + Poc div v - Q (U) 


R 


M, 


oo 


P T 


( 12 ) 


Equations (7) to (12) are the equations of conservation of mass, momentum and 
energy, respectively. Equation (12) is the equation of state. It is 


(9, y, r) spherical coordinates 
t time 

v = (u, v, w) vector of wind velocity with Es components in south-, east- 
and upward direction, 
p pressure (p 00 mean pressure) 

P density (p 00 mean density) 

T temperature (T 00 mean temperature) 

0. angular frequency of the earth's rotation 
v number of collisions between the ions and one neutral 
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I Inclination angle of the earth's magnetic dipole field 

g magnitude of gravitational force 

c v specific heat at constant volume 

k 00 mean coefficient of heat conductivity 

R gas constant 

M 00 mean molecular weight 


div v - 


1 3 1 

7^030 < sin£? u ) + r sin 0 3IV 


3 v 1 3 „ 

+ (r 2 w) 

r 2 3 r v ' 


A = 


1 3 2 1 

-f 

r 2 sin 2 G 3^ 2 r 2 sin Q 


Te ( sind Te) 


4. Development of the Physical Parameters in Series of Spherical Harmonics 

Our model shall be valid during equinox conditions. Therefore all asym- 
metric spherical functions in the development of density, temperature and pres- 
sure disappear. Because of our assumption in Equation (6) higher order terms 
of the coefficients in the series of spherical harmonics are neglected. We use 
the unnormalized associated Legendre polynoms, which are 

P oo " 1 

P 10 = cos Q ; Pjj = sin 6 

P 20 - ~ 2 ( 3 cos 2 9 - 1) ; P 21 = 3 sin 9 cos 9 

We start with the expressions for temperature, density, pressure and vertical 
wind. They have the general form 
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a(r, 9, \) 


(13) 


a oo( r ) p oo + a 20 ( r ) P 20 (^) + a n ( r ) p i i (^) e ' 


j Qr 


Here, of course, the mean vertical wind disappears (w 00 = 0). Because a n (r) 
is complex, we consider the real part 


2 Real j a u (r) 
as the physical solution. 


- |A U | cos 


n ( r “ r u) 


|A U | = 2|a n | 


is the magnitude and 


is the time lag. 


li 


ar g(a u ) 

n 


It is 

\ 

r = t+ lT 

the local time, while t is the universal time. 

In the development of horizontal winds into series of spherical harmonics it 
is convenient to add a factor (sin 9) ±l to the series. This allows a direct trans- 
formation of the spherical harmonics of equal domain m within the system of 
Equations (7) to (12), merely by using the well known recurrence formulae of 
spherical functions. 

We have to introduce a mean longitudinal wind v 00 in order to obtain a 
unique solution of the problem. Meridional winds as well as longitudinal winds 
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are symmetrical with respect to the equator during equinox. Longitudinal winds 
and the zonal-component of the meridional wind disappear at the poles. But the 
time dependent component of the meridional wind is different from zero there, 
because a finite pressure gradient can be maintained at the poles. 


In order to fulfill the conditions outlined above we introduce the expressions 


u(r, 8, k) - sin0u lo (r)P lo (0) + P 21 (0)e^ T 


(14) 


v(r ,6,k) - sin0{v oo (r)P oo + v 20 (r)P 20 (6') + v 11 (r)P u (0)e jfiT } (15) 


We have reasons to expect that this representation of the horizontal winds 
possesses the fastest convergence rate, thus allowing as a first order approxi- 
mation the negligence of higher order terms in the series. 

We furthermore assume known functions of altitude of the coefficient of heat 
conductivity /< 00 = A(r) yT 00 (r) and of molecular height M 00 (r). Likewise, the 
EUV heat source Q is assumed to be developed in a general form of Equation (13). 
As one can easily see from a harmonic analysis e.g. of a Chapman-function of 
the ion production the three coefficients in Equation (11) constitute the predomi- 
nant heating terms. 

Finally, we take a known function of the collision number 

V ~ V 00 ^ r ) ^*00 + V 20 ( r ) ^20 (^) ( 1 ®) 

because the coefficients v 00 and v 20 have comparable magnitudes. The angle of 
incidence I of the geomagnetic dipole is related to the colatitude in a well known 
manner. We obtain 
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cos 2 I 


sin 2 9 


sin 2 I > = < 4 cos 2 9 


sin I cos I 


2 sin 9 cos 9 


(1+3 cos 2 0 ) 


and develop 


“oo^no + a on^ > on^^ + 


1 + 3 cos 2 6 00 00 20 20 


With 


a oo 


0.604 


(Xjo -0.520, 


neglecting higher order terms. 


With these assumptions we enter Equations (7) to (12), use the recurrence 
formulae of the spherical harmonics (see e.g. JahrJke and Emde, 1945), collect 
all coefficients belonging to the domain (n, m) and ob' 'in for each Equation (i), 
numbered according to the Equation-number (7) to (1. ; expressions of the general 
form 

CO 00 

YL F ™'' *' 0 ( 17 > 

n = 0 m = 

Because of the orthogonality of the spherical harmonics each coefficient F< * > 
must be zero. Thus, we find the set of equations (for abbreviation we write: 
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The right hand side of Equation (18a) is not exactly zero but consists of 
terms of the same order of magnitude as the terms in Equation (19f). As one 
can easily see, these terms are two orders of magnitudes smaller than the terms 
on the left hand side in Equation (18a). If they are neglected, Equation (18a) is 
the barometric height formula. Now the system of Equations (18) is completely 
decoupled from the system of Equations (19) and gives the mean values averaged 
over time and sphere. 

Equations (19) belonging tc the domain m = 0 describe the time independent 
latitudinal variations of the thermospheric parameters. This system will be 
called the zonal system. Equations (20) belonging to the domain ir. = 1 give the 
diurnal variations. We call this domain the diurnal system. 

5. Gener a l Solution 

The approximate validity of the barometric height formula [Equation (18a)] 
allows the determination of mean temperature, pressure and molecular weight 
from the observed mean density by a static diffusion model (Jacehia, 1964). In 
our treatment of the diurnal tide within the thermosphere we therefore staid; with 
these known parameters of the thermosphere and use them for the normalization 
of the set of Equations (19) and (20). We write for abbreviation 
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w<°> = w 20 /C 

p(°) - P 20 /p 0 o 

T(0) = t 20 /t 00 

y (0) = ^oono/^Poo) 

V (0) = v 00 /C 
v«» = v 2Q /C 

u (0) = U i 0 /C 

p^°) - /0 20 /p 00 

Q (0) = Q 20 /(»Poo) 


w (l > = w u /C 
P (1) = Pi l' Poo 
T(1) = T u/ T oo 

y (1) _ K oo 1*11/ (^Poo ) 

- Pn/Poo 
v< x > = v n /c 
u (x ' = u 21 /C 

Q (1) = Qu/(fiPoo) 


( 21 ) 


We eliminate the parameters p (i ^ , v (i) , u (i) and V (0) from Equations (19d) 
to (19g) and (20d) to (20f), respectively, and obtain two independent systems of 
first order linear ordinary differential equations, which we write in concise 
matrix form 


d»<*> 


= K< x >(£) e< x > + h^>(£) 


( 22 ) 
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with 
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(i) 


p(i ) 


h<i) 


) 



r Q(i) 


; (i = 0,1) 


and 


£ 


n 

c 


r 


C = velocity of sound 

K( 0) (g) is a 4 x 4-matrix with real elements, K (1) g) is a4 x 4-matrix with 
complex elements. The elements of K (i ) can be found from Equations (It/) and 
(20), respectively. The elements of K( ; ) as well as of e. are dimensionless 
and are of comparable magnitudes. 

Equation (22) can be solved by standard methods. If we approximate the 
real atmosphere by a number of isothermal homogeneous slabs of thickness A£ r 
in which the elements of K <* > are constant, then the solution of Equation (22) 
is 


eGf n ) = P 0 ne ^o) +<t> 0 


(23) 
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The matrix r. matches the boundary conditions of continuous wave parameters 

w, p, T and k T* at the boundary between two adjacent slabs of temperature T. 

and T. . 

1 


22 



6. The Horizontal Wind System 


From Equations (19e) to (19g) and from (20e) and (20f), respectively, we 
notice that the horizontal winds are connected with the vertical wind and the 
pressure by linear equations. Since the vertical wind is at least one order of 
magnitude smaller than the horizontal wind components, we can neglect the 
vertical wind. In the following discussion we shall neglect furthermore the 
higher order terms v 2Q and a 20 of the collision number and of the magnetic dip 
angle, because they influence the numerical values of the winds but not their 
general behavior. In the numerical calculations which will follow in the next 
section we shall however take into account these higher order terms. 


Negligence of vertical winds and of latitudinal dependence of collision num- 
ber and dip angle leads then to very simple relations between the horizontal wind 
components and the pressure field. These are for the zonal system (m = 0) 
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For the diurnal system (m = 1) we obtain 
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For large collision numbers (v 00 /fi >> 1) it is 
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For small collision numbers (v 00 /f2 << 1) only the meridional zonal velocity 
u io a v oo depends on v Q0 . The magnitude of the horizontal winds therefore is 
very sensitive with respect to the collisions between neutrals and ions . 
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The longitudinal wind of system m = 0 is 


v zonai " sin0 ( v oo P oo + v 2 o P 2o) “ 3 sin 6 cos 2 6 V Q() (26) 

This wind is the geostrophic wind due to the meridional pressure gradient p 20 . 
Collisions between neutral and ions cause a mean meridional wind u 10 to flow 
in the direction of the pressure gradient. Collisions moreover reduce the mag- 
nitude of the geostrophic wind v zonal . 

A mean longitudinal wind [Equation (26)] means that the whole atmosphere 
rotates with respect to the earth surface. It has been pointed out by King-Hele 
and Allan (1966) that such wind exists within thermospheric heights. King-Hele's 
winds, derived from the observation of the changes in the orbit inclinations of 
satellites, blow from west to east and are of the order of 100 m/sec. One of 
King-Hele’s explanations for the origin of these winds is in fact the geostrophic 
wind hypothesis. The pressure gradient built up to generate this wind must be 
directed toward the poles. Observation of Jacchia and Slowey (1967) as well as 
density measurements on board of Explorer 32 (Newton, 1968) however 
show that during low solar activity (F < 150) the zonal pressure gradient is 
directed toward the equator above 300 km altitude (p 20 > 0). According to 
Equation (24) this pressure gradient causes a meridional wind u 10 blowing 
toward the equator which is in agreement with a hypothesis to explain the equa- 
torial F2 anomaly (Mayr and Volland, 1969). At moderate solar activity (F - 150) 
the zonal pressure p 20 disappears and it seems likely that it becomes negative 
at high solar activity. According to Equations (24) and (26) we therefore expect 
a zonal wind blowing toward the west at low solar activity and blowing toward 
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the east at high solar activity. Gexsler (1967) too finds this westward zonal wind 
at solar minimum conditions from his numerical study of the horizontal wind- 
system in the thermosphere. King-Hele however states from his observations 
that the mean wind always blows toward the east regardless of the solar activity 
with even a greater magnitude at solar maximum than at solar minimum. 

There exist several possibilities to explain this discrepancy between King- 
Hele’s observations and the theoretical results. First, it will be shown in 
Figure 2 that due to the higher order terms of the collision number and the geo- 
magnetic dip angle the calculated zonal wind blows toward the east within ± 15° 
from the equator at 300 km altitude. Second, we neglected in our calculations the 
time dependence of the collision number. A diurnal term v ix would couple the 
zonal system with the diurnal system. It has been pointed out by Voiland (1966), 
that such coupling would in fact give rise to an additional eastward component of 
the mean wind. Third, electric fields, entirely neglected in our treatment, could 
cause ion drift such that the observed mean wind results (Hines, 1965). 

We should bear in mind however that the zonal pressure gradient not only 
depends on solar activity but also on height. For a better understanding of this 
problem we need numerical calculations of the three dimensional wind and pres- 
sure fields, and furthermore we have to wait for more detailed observational 
data of the zonal wind depending on height, latitude and solar activity before a 
final answer can be given. 
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7. Numerical Cal culations of the Horizontal Wind Field 

For a numerical study of the horizontal wind field at 300 km altitude we used 
the following data valid at low solar activity (F = 100): 


Mean temperature and molecular weight from CIRA model 3. From (Volland 
et al., 1969) the relative diurnal pressure 
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From Newton (1968) the zonal pressure, extrapolated from his density and 
temperature data 

P 20 

— = 0.25. (28) 

p oo 

The collision number from the maps of the F2 critical frequency (Martyn, 1955) 


v 00 = 2.4 x 10 -4 sec -1 

(29) 

v 20 - -1.9 x 10" 4 sec -1 . 


We calculated the horizontal winds from Equations (19e) to (19g) and (20e) and 
(20f), respectively, and neglected the vertical wind components. Figures 1 and 2 
give the calculated zonal winds 


u zona i = sin9 u io p io( 0 ) 

V zonal = sin0|v oo P oo + V 2O P 2O (0)l, 


(30) 


27 




the mean longitudinal wind 
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and the magnitudes of the diurnal winds 
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versus latitude. Moreover, the total horizontal winds 
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where the time lags have been calculated to 


r j 1 = 21 00 local time 
T = 2 20 local time, 

are plotted for the time of maximum and minimum winds respectively. 

We notice that due to the higher order terms v 20 and a 20 of the collision 
number and of the geomagnetic dip angle the zonal longitudinal wind v zonal has 
a small eastward component near the equator (Figure 2) . 

Furthermore, the meridional wind u blows at low latitudes toward the 
equator even at day time. This unexpected result- follows from the relatively 
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large zonal u t0 -component which predominates the meridional velocity at low 
latitudes. Such equatorial wind at noon can in fact create the equatorial F2 
anomaly: the relative maximum of the F2 electron density at ±10° latitude and 
the trough at the equator. The electrons drift upwards by the equatorial wind 
and cause an enhancement of the electron density near the equator (Mayr and 
Volland, 1969). 

The wind u 10 depends on the zonal pressure gradient p 2Q . ^his pressure 
gradient disappears at heights above Aoc-fom at moderate solar activity and 
probably becomes negative with further increasing solar activity (Newton, 

1968). As a consequence the equatorward directed wind at noon decreases or 
even reverses with increasing solar activity, giving rise to a weakening or disap- 
pearance of the F2 anomaly. This picture is in qualitative agreement with the 
observations. 

An equatorward wind implies of course a total pressure gradient directed 
toward the equator. The combination of the pressure terms p 20 and p n chosen 
in Equations (27) and (28) therefore gives rise to a small relative pressure mini- 
mum at the equator at noon. The elongation of the density bulge at noon and at low 
solar activity observed by Jacchia and Slowey (1967) points in fact into this direc- 
tion. It is possible that due to the low degree of space and time resolution in the 
satellite drag measurements this minimum has been overlooked. 

Whether the p 20 term is mainly created by solar EUV or by another heating 
mechanism (e.g,: the regular component of the precipitation of high energy 
electrons into the auroral zones related to the geomagnetic S D -current) is an 
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open question. The different components of the pressure gradient are certainly 
functions of height. Since most of the heat is deposited within a height range 
below 200 km we expect a pressure gradient built up there and directed away 
from the heating zone. In the case of solar EUV heating this would be a gradient 
toward the poles. Due to this pressure gradient a global wind system is set up,, 
in which the vertical wind plays a decisive role. There must be a height range 
— probably the region above the main heating zone — where a return flow achieves 
flow continuity. This return flow however must be connected with a pressure 
gradient directed opposite to th pressure gradient within the main heating zone. 
Because of the important role of the vertical wind a quantitative treatment of 
such problem is only possible by a numerical integration of the dynamic equa- 
tions of the spherical thermosphere [Equations (23)] . An attempt to determine 
the vertical winds from the pressure field of the Jacchia model and from the 
horizontal winds of Geisler (1967) via the equation of continuity (Dickinson and 
Geisler, 1968) suffers from the unaccuracy of the density observations below 
200 km altitude. 

Figure 3 finally presents the horizontal wind field in the northern hemisphere 
at 300 km height in a form used by Kohl and King (1967). This allows a direct 
comparison between Kohl and King's winds (their Figure 5) and our wind. The 
striking difference between both wind fields is the presence of the equatorward 
wind component at low latitudes and of daytime in our wind field. This component 
is due to the relatively large zonal wind. Its generating zonal pressure gradient 
p 20 can be observed already in Jacchia and Slowey's (1967) data at Jow solar 
activity, though Jacchia's number value is smaller than the value observed by 
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Newton (1968). Kohl and King (1967) entirely neglected the zonal winds. They 
moreover used a constant collision number for the whole hemisphere. 

Therefore we can only compare our diurnal windsystem (u 21 , v n ) with their 
data, which is in reasonable agreement. It shows that slight differences in the 
numerical data of the pressure gradient p u and of the collision number v are 
relatively unimportant compared with the influence of the zonal winds. 

Geisler's (1967) winds are very similar to Kohl and King's (19671 winds 
because they are based on the same calculation method. Gei iC-r, however, found 
already the mean component v 00 from Jacchia's data. Since Jaoo.._ 's zonal 
pressure component is supposed to be too small, Geisler could not detect the 
equatorial winds at daytime, though he observes "a convergence in low latitudes" 
in his horizontal wind field, which is the direct consequence of the superposition 
of the diurnal system with the zonal system. 

8. Conclusion 

The solar diurnal tides of the spherical thermosphere have been treated 
using spherical harmonics of lowest degree as eigenfunctions of the problem. 
Perturbation theory is a sufficient approximation within the thermosphere below 
400 km altitude. It allows to determine the mean physical parameters like pres- 
sure, temperature and molecular weight from the observed density data by a 
static diffusion model (Jacchia, 1964). Thus, these data were the basis for a 
perturbation treatment of the hydrodynamic equations leading to a complete 
separation between a zonal component, depending only on latitude and height, and 
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a diurnal component, dependent on latitude, height and local time. A general 
solution for both the zonal and the diurnal system is given. 

The horizontal wind field turns out to depend mainly on the horizontal pres- 
sure gradients. Based on observed data for the zonal and the diurnal components 
of the pressure gradient, the horizontal wind at 300 1cm altitude and for low solar 
activity has been calculated. Its results show an equatorward directed wind flow 
at low latitudes even at day time which is due to the relatively large zonal pres- 
sure gradient. This wind could possibly cause the equatorial F2 anomaly: the 
trough in the electron density at the geomagnetic equator during low solar 
activity. 
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Table I 

Number S nm of the _i gen functions of the diurnal tidal component 

(m = 1) for two different isothermal atmospheres of temperature 
T 0 = 250°K and T 0 = 1000°K. 



P n 

P 3l 

P 5 1 

T 0 - 250°K 

0.88 

2.16 

3.41 

T 0 = 1000°K 

2.16 

' 

5.29 

8.36 
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LATITUDE 
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figure 1. Components of calculated latitudinal wind u 
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